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Abstract. 

We consider energy deposition of high energy electrons and photons in universe. 
■ We carry out detailed calculations of fractions of the initial energy of the injected 

electron or photon which are used to heat, ionize and excite background plasma in the 
^ ■ early universe for various ionization states and redshifts. 

o 

c/3 , 1. Introduction 

The energy deposition of fast electrons and high energy photons in partly ionized 
^ ■ plasma is an important issue in the wide range of physics and astrophysics. In 
cosmology, for example, the high energy photons and/or electrons are injected from 



X 



G^ ' decay or annihilation of massive particles. The recent observation of WMAP on Cosmic 
^ ■ Microwave Background (CMB) has shown that dark matter constitutes about 80% of 
■ the total mass in the universe [3]. One of the most promising candidates for dark 
^ ' matter is a supersymmetric particle with mass O(IOO) GeV. If such massive particles 
are annhilate and/or decay, the ejected charged particles and photons interact with 
various background particles: electron, atomic hydrogen and CMB photon, and hence 
gives the significant effects on the thermal histrory of the universe. In this case, to 
estimate the cosmological effect of dark matter annihilation/decay precisely, we need 
to understand development of electromagnetic cascade showers induced by primary 
particles, energy loss of charged particles and photons, and so on. Besides dark matter, 
there are many candidates for inducing electromagnetic showers and affect the cosmic 
background plasma in physics beyond the standard model. 

Detailed calculations for electrons and photons slowing down in partially ionized 
plasma of atomic hydrogen were carried out by many authors (see references in [1]). 
However, these authors mainly studied the case of injected electron energy up to keV. 
In this paper, we extended these calculations which are valid up to much higher energy, 
i.e., 1 TeV. The injected high energy particles produce the cascade showers and lose their 
energy. We can categorize three types of energy loss: heat, excitation and ionization 
according to what the energy of the particle is used for. It is the aim of this paper to 
derive these quantities precisely. They depend on the energy of primary particles as well 
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as redshift and state of ionization. There are two reasons why it is difficult to derive 
them. One is the large number of interactions involved with the energy degradation. 
The other is connection between the energy degradation of charged particles and that 
of photons. For example, high energy electrons can produce photons as much the same 
energy through inverse Compton scatterings with CMB photons. In other words, we 
can not calculate the evolution of the energy of charged particles and photons separately 
but simultaneously. 

In section [2], we show the numerical method to calculate the energy degradation 
of the primary particles. In section [3l the evolutions of the energy of primary electron, 
photon and positron are presented. The results of the energy degradation of the primary 
particles is shown in section |H In section [5], we summarize the results. 

2. Numerical Method 

For incident electrons and photons, there are many interactions which contribute to 
energy degradation. It is convenient to divide these processes into two groups. One is 
the group characterized by losing only a very small portion of energy in one collision 
(continuous loss). The other is the group characterized by being likely to lose a 
significant portion of energy in one collision ('catastrophic' loss [6]). In the case of 
the latter, it is necessary to use an integro-differential equation to calculate electron 
energy spectrum. 

Let El, E2, ■ ■ ■ , E]y be a discrete set of energies of particles {Ei < Ej for i < j) and 
Np{Es)AEs be the number of particles with energy between Es—AEs/2 and Es + AEs/2. 
The accuracy of the numerical method is limited by the bin size (AEi). Since we consider 
very large energy region ranges from 10 eV to 100 GeV, the bin sizes are taken to increase 
as energy so that AEi/ Ei become constant. The particle energy spectrum is given by 



where P{Ei, Ej)dt is the probability that a particle with energy Ei undergo a collision 
causing it to lose energy and have energy Ej in time dt, and Q{Ei) represents sources and 
sinks of particles corresponding to possible production, annihilation or gradual leakage 
from the energy range that we consider [1]. The ffist term of the R.H.S. of Eq. ([T]) 
represents continuous loss and the second and third terms represent outflow and inflow 
caused by catastrophic loss. 

In this paper, we are mainly interested in how much initial particle energy convert 
to heat, excitation and ionization. For convenience, we define Xh{E), Xe{E) and Xi{E) as 
fractions of the initial energy E which go to heat, excitation and ionization, repectivelyi^] 

I Wc frequently use "fraction of heat" , "fraction of excitation" and "fraction of ionization" to refer Xh , 
Xe and Xij respectively. Please do not confuse Xi with "ionization fraction" Xe which is the fraction of 
ionized hydrogens. 





(1) 
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Suppose that we know already Xh{Ei), Xe{Ei) and Xii^i) with energy less than 
Eg and consider a particle with initial energy Eg. Since there are no particle with 
energy more than Eg, the energy degradation is characterized by the first term of the 
R.H.S. of Eq. ([1]) (continuous loss), the second term (catastrophic loss) and the last 
term (annihilation). Here we only consider that annihilation of the primary particle 
contributes to Q{Es). All these terms are in proportion to Np{Es), so we define the 
following number loss function 

' -dE' 



L{E, 



dt 



+ J2ntVpai{Es) (2) 



where Ut is the number density of target particle, Vp is the particle velocity and <Ji{Es) 
is the catastrophic and annihilation cross sections. From number loss function, we can 
obtain the probability that a particle undergoes a particular collision "m" [5]. The 
collision frequency for a particular continuous loss may be defined by 

-dE 



AE. 



dt 



(3) 



and the collision frequency for a particular catastrophic loss "a" and annihilation 
are given by 

UaiEs, Ei) = UtVpa^iEs, Ei), (4) 
up{Es) = ntVpap{Es). (5) 
The total collision frequency for the catastrophic loss is given by 

Ua{Es) = I dEntVpaa{Es,E) (6) 



Then the probability P{Es, Ei) is written as 

Em f^miEs, Es-l) + 22aMEs) + MEs)) 
En.'^m{Es.E,_,)+ZMEs) + MEs)y 

Combining these probabilities with the data about Xh{E), Xe{E) and Xi{E) with energy 
less than Es, Xhi.Es), Xe{Es) and Xi{Es) can be obtained. Please notice that the 
definition of the frequency for continuous loss depends on the size of the energy bin 
but final result is independent of it. The reason is as follows. When the all bins are 
divided into halves, the probability P{Es, Ei) decreases by fifty percent if the bin size is 
small enough. This is because the collision frequency for continuous loss is much larger 
than that for catastrophic loss in this case. However, the effect of the decrease of the 
probability is cancelled by the increase of the bin numbers. We have checked that our 
results are independent of the size of the energy bin. 
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Figure 1. Continuous energy loss rates for electrons with Xg ~ 10~^ and 1 + z = 1000. 
Solid line represents elastic collision with atomic hydrogen, dashed line represents 
Coulomb loss with background electron and dotted line represents inverse Compton 
scattering with CMB photon. 
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Figure 2. Cross sections for electron impact excitation and ionization of H. 



3. Individual Evolution 

3.1. Electron 

For incident high energy electrons, there are many interactions which contribute to 
energy degradation: elastic collision with atomic hydrogen, excitation and ionization 
of atomic hydrogen. Coulomb loss with background electrons and inverse Compton 
scattering with background photons (CMB photons). The cross sections and energy 
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Figure 3. Electron energy degradation with = and 1 + z = 1000. Solid line 
represents Xh, dashed line represents Xex and dotted line represents Xi- Dot dashed 
line represents the fraction of photon energy which energy is larger than Ry. 



loss rates of these interactions are described in Appendix A and shown Figs. [T] and O 
Hereafter, we use electron kinetic energy (Kg) instead of electron energy [Ee = Ke + rrie) 
for convenience. In this paper, we have used 4 level (2s, 2p, n = 3 and n > 3) 
approximation in considering electron impact excitation. From Fig. [1] it is found that 
energy loss is dominated by collisions with background electron (Coulomb loss) at low 
energy and inverse Compton scattering off CMB photon at high energy. The reason is 
simple. At low energy, the average energy loss of electron is roughly AE^ oc P'^Ecmb 
for one inverse Compton scattering. So electrons lose very few fraction of their energy 
and hence energy loss is dominated by collisions with atomic hydrogen and background 
electron. On the other hand, at high energy, since the number density of CMB photon 
is much larger than that of background electron, energy loss is dominated by inverse 
Compton scattering. Energy loss by Coulomb collisions is so efficient that free electrons 
with ionization fraction, Xg, exceeding 10~^ have a substantial influence on the energy 
degradation [1]. Thus, as the ionization fraction increases, the fraction of the initial 
electron energy which is converted to heat {xh{E)) increases as shown later. 

We calculate the energy degradation of electron following the method described in 
the previous section. There are two free parameters: ionization fraction Xe and redshift 
1 + z. Ionization fraction is relevant for collisions with atomic hydrogens and background 
electrons, and redshift is mostly relevant for inverse Compton scattering. For simplicity, 
we regard the distribution of CMB photons as monoenergetic {E^ = 6.34 x 10~^'^(1 + z) 
GeV), not blackbody. Free and bound proton number densities can be parameterized 
by baryon-to-photon ratio rj: Up = rjUcMB- We adopt rj = 6.1 x 10"^*^ from the result 
of WMAP 3-year observation [3]. We take into account properly secondary electrons 
which is produced by electron impact ionization. 
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Figure 4. Electron energy degradation with ~ 10^^ and 1 + z = 1000. Solid line 
represents Xh, dashed line represents Xex and dotted line represents Xi- Dot dashed 
line represents the fraction of photon energy which energy is larger than Ry. 




Figure 5. Electron energy degradation with Xe ~ and 1 + z = 100. Solid line 
represents Xh, dashed line represents Xex a-nd dotted line represents Xi- Dot dashed 
line represents the fraction of photon energy which energy is larger than Ry. 



Energy degradation for several values of Xe and 1 + 2 is shown in Figs. [3H5l Here, 
we plot Xh{E), Xex{E) and Xi{E). In addition, we also plot the fraction of the initial 
energy that goes to photons with energy larger than Ry(= 13.6 eV). This energy is 
the threshold energy for photoionization, and it is granted that these photons also 
contribute to heat, excitation and ionization and we will treat them in Section HI On 
the other hnad, photons with energy less than Ry only heat background electrons with 
Compton scattering. An oscillating behavior below 50 eV in Figs. [3] and [5] reflects 
the nature of discrete energy loss. This behavior is not seen in Fig. H] since the effect 
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of Coulomb loss dominates over ionization and excitation. The common features in 
Figs. [3ll5] are as follows: (1) the fractions of heat, excitation and ionization are the 
same order for Ee < 10^ eV, (2) heat dominates at IC^ < E^ < 10^ eV, and (3) finally 
the energy of photons dominates for E^ > 10^ eV. At relatively low energy, inverse 
Compton scattering is inefficient. Moreover, the cross sections and energy loss rates 
for momentum loss, exitation and ionization are almost the same order if ionization 
fraction is not large, i.e. Coulomb loss is not dominant. As a consequence, Xh, Xex 
and Xi ci-re nearly the same order and our calculation corresponds to the result of [1] 
up to Ef. ~ keV. As electron energy increases, inverse Compton scattering becomes 
significant. If the energy of scattered-up photon is less than the threshold energy of 
photoionization, the energy loss due to inverse Compton scattering converts to heat and 
the fraction of heat approaches to unity. When electron energy further increases and 
most of scattered-up photons have enough energy to ionize atomic hydrogen, the energy 
of incident electron exclusively converts to the photon energy. 

Next we will show how these fractions are dependent on parameters. To see the 
effect of ionization fraction, let us compare Fig. [3] with Fig. HI It can be seen that the 
fraction of heat increases at low energy as ionization fraction increases. This is because 
Coulomb collision converts initial electron energy into heat exclusively. Besides, the 
collision frequencies for ionization and excitation are proportional to the number density 
of hydrogen atom, (1 — Xe), and hence they are alomost independent of if Xg <^ 1. 
Thus, Xex{E) and Xi(-E) is smaller as ionization fraction increases. At high energy, the 
results are irrelevant for ionization fraction since inverse Compton scattering dominates 
over all other processes. Comparing Fig. [3] with Fig. [5], it can be seen that the rise of 
heat and photon energy shifts to high energy side as redshift decreases. This is because 
the effect of inverse Compton increases in proportion to the energy of CMB photon, and 
hence in inverse proportion to redshift. 



3.2. Photon 

As well as electron, there are many interactions which contribute to energy degradation 
for incident high energy photons: photoionization, Compton scattering with background 
electrons, pair production in matter, photon-photon scattering and double photon pair 
creation. Here, we neglect photo-excitation since the resultant excited state immediately 
emits photon and goes down to the ground state. The cross sections and energy loss 



rates of these interactions are found in Appendix A and shown in Fig. [61 



Energy degradation is plotted for Xg = and Xe = 0.98 in Fig. [7] and Fig. [HI 
Here, we plot Xh and Xi- addition, we also plot the fraction of initial energy which 
goes to electrons with energy larger than 0.75Ry. This energy is the threshold for 
electron impact excitation of atomic hydrogen. Electrons with lower energy only heat 
background electrons through elastic collision with atomic hydrogen and Coulomb loss. 
The sudden falls of Xh reflects the nature of discrete energy loss. The energy for the flrst 
fall {E.y = Ry) corresponds to the threshold energy of photoionization and the second 
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Figure 6. Photon collision frequencies with ~ and 1 + z = 1000. 
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Figure 7. Photon energy degradation with Xe = and 1 + z = 1000. Solid line 
represents Xh, dashed line represents Xi ^-nd dotted line represents the fraction of 
electron energy which energy is larger than 0.75Ry. 



one {Ey = 7/4Ry) corresponds to the sum of the threshold energy of photoionization 
and that of electron impact excitation. In that case the energy of the incident photon 
goes into ionization and the electron produced in the ionization, and is not used for 
heating. (Of course the electron further interacts with background plasma and heats it. 
This will be taken account in the next section.) The rise of Xh around E^ ~ lO^eV is 
due to Compton scattering. Unlike the case of incident electron, the fraction of electron 
energy dominates even at low energy. This is because energetic electrons are produced 
through photoionization at low energy and Compton scattering at high energy. 

The increase of ionization fraction causes the increases Xh especially at low energy. 
This is because photoionization and Compton scattering determine the degradation 
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Figure 8. Photon energy degradation with 1 — = 10~^ and 1 + z = 1000. Solid 
line represents Xh, dashed Une represents Xi ^.nd dotted hne represents the fraction of 
electron energy which energy is larger than 0.75Ry. 



of electron energy at low energy. The effect of photoionization decreases as ionization 
fraction increases, and hence Xh increases. Since the photoionization rate is proportional 
to (1 — Xe), the effect of photoionization becomes small as Xe approaches 1, which leads 
to relative enhance of Compton scattering and increases Xh- Note that we use the the 
baryon density instead of the electron density when we calculate the energy losses due to 
Compton scattering. This is because Compton scattering only becomes important when 
is sufficiently larger than Ry and the interaction is insensitive to whether an electron 
is bound or not As redshift increases, photon-photon scattering and double photon 
pair creation become important. These processes, however, become dominant at very 
high energy as seen in Fig. [61 This effect is not appeared until the energy deposition of 
electron and positron taken into account. 

3.3. Positron 

Positrons are produced by pair production in electric field of nuclei and double photon 
pair creation. Therefore, we should take the energy degradation of positrons into 
account. Although the energy degradation of positrons is almost the same as that of 
electron, there are two differences between electron and positron. One is the sign of its 
charge. The other is the indistinguishability between a primary electron and a secondary 
electron in the processes where target is an electron or an atomic hydrogen. These 
differences become less important at high energy. For simplicity, we assume positrons 
lose their energy just like electrons in this paper. Besides, positrons finally annihilate 
with background electrons through either free annihilation or the formation and decay of 
positronium [30]. It depends on the temperature, density and state of ionization of the 
background electron which process is dominant. Roughly speaking, most of the positrons 
undergo annihilation after the significant loss of its energy. Therefore, we assume that a 
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Figure 10. Photon energy degradation witli Xe = 10 and 1 + z = 1000. Solid line 
represents Xh- dashed line represents Xi and dotted line represents Xex- 



positrons forms a positronium with a background electron after losing almost all of its 
energy and decay. The positronium annihilates to two photons (each 0.511 MeV) 25% 
of the time, and to three photons (each less than 0.511 MeV) 75% of the time. The 



energy spectrum from the three-photon annihilation is described in Appendix A 
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4. Results 

As previously mentioned, we should calculate the evolution of the energy of charged 
particles and photons simultaneously. In Fig. [9] and Fig. [10], we show the fraction 
of heat, excitation and ionization when the primary particle is electron and photon 
respectively. Let us examine these figures. 

First, we consider the electron case. When the energy of the electron is small 
[Ke < 10^ eV) and the ionization frcraction is not large, the fractions of heat, excitation 
and ionization are roughly the same order. This is because relative smallness of 
ionization fraction makes the effect of Coulomb loss comparable with that of excitation 
and ionization of atomic hydrogen. The fraction of heat increases as ionization fraction 
increases and vice versa. For larger electron energy {K^. > 10^ eV) , the dominant energy 
loss mechanism is the inverse Compton scattering. If the electron is non-relativistic, 
the energy of scattered-up photons is so small that these photons only contribute to 
heating of the background particles. As a consequence, the fraction of heat reaches 
near unity. When an electron is relativistic, the energy of scattered-up photons exceeds 
the threshold energy of ionization of atomic hydrogen. In this case, it is necessary to 
estimate the contribution of these photons. It is seen that the fractions of excitation 
and ionization are a little larger than that of heat with low photon energy in Fig. [101 
For this reason, all of these fractions becomes almost the same amount. The oscillating 
structure of the x's around — 10^ ~ 10^° eV reflects that of scattered-up photon 
around ^ 10^ ~ 10^ eV. When the energy of electron becomes ultra- relativistic, 
these oscillating structure is averaged out and vanishes. 

Next, let us consider the photon case. When the energy of the incident photon 
is small {E^ < 10^ eV), the fractions of excitation and ionization are larger than that 
of heat. This is because the dominant energy loss mechanism is photoelectron effect 
unless ionization fraction is very close to unity. In this case, a photon ionizes an atomic 
hydrogen and emits a photoelectron whose energy is almost the same as the primary 
photon. Therefore, the behavior of photons is very similar to that of electrons at low 
energy. As photon energy increases, Compton scattering becomes the dominant energy 
loss. When the photon energy is larger than the electron mass, the forward scattering 
becomes dominant and the energy of the recoil electron approaches to that of the photon. 
In other words, Compton scattering produces a recoil electron with energy same as the 
incident photon. Therefore, the behavior of photons traces that of electrons in this 
energy region. 

When photon energy is larger than 10^ eV, the effect of pair-production in 
matter, photon-photon scattering and double photon pair creation can not be neglected. 
However, it is hard to explain the influence of these processes because the most dominant 
process among the three depends on redshift and photon energy, and the shape of the 
spectrum of secondary particles is very sensitive to primary photon energy. At high 
energy, however, the fraction of heat, excitation and ionization become constant values 
just like electron. As high energy photons are closely related with electrons, these 
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Figure 11. Electron and photon energy degradation with 1 + z = 1000. Sohd hne 
represents Xh, dashed hne represents Xi f-nd dotted hne represents Xex- 
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Figure 12. Electron and photon energy degradation with Xe = 10 Solid line 
represents Xh, dashed line represents Xi and dotted line represents Xex- 



constants will become the almost same value in either case. 

To see how ionization fraction Xe and redshift z have influence on the energy 
degradation, we plot the fractions of heat, ionization and excitation as function of 1 — 
and 1 + z in Figs. [H] and [121 repectively. Here we have taken 10^^ eV for the initial 
energy of electrons or photons. This is because all the fraction become constant values at 
this energy. Unless ionization fraction is very small, Xh and Xi are almost independent 
of it. This is because Xi is mainly determined by photoionization due to low energy 
secondary photons in this case and the collision frequency of photoionization is much 
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larger than that of Compton scattering at low energy [Fig. E]. Therefore the effect of the 
change of ionization fraction is almost irrelevant. However, the fraction of excitation Xe 
is very sensitive to ionization fraction since the collision frequency of electron impact 
excitation is in proportion to 1 — Xe while the competing processes (inverse Compton 
and Coulomb loss) are independent of 1 — Xg. The fraction of heat Xh does not depend 
on 1 — since Xi is independent of 1 — Xe and Xex is very small quantity. (Please notice 
that Xi~^Xex + Xh = !•) When the ionization fraction is very small, the effect of Coulomb 
scattering is weakened and these fractions become the same order. Unlike ionization 
fraction, the change of redshift seems to have no influence on the energy degradation. 
This is because the change of redshift effects which process is dominant at high energy, 
but is almost irrelevant at low energy. Xi ^iid Xe are determined by secondary low energy 
particles and these fractions are nearly independent of redshift. 

5. Conclusions 

We have carried out detailed calculations of the fractions of the initial energy of the 
injected electron or photon which are used to heat, ionize and excite background plasma 
in the early universe. In the high energy limit [E > 10^^ GeV), we have shown that the 
fractions are alomost independent of the initial energy. Our calculations are valid up 
to TeV and can be applied to various cosmological and astrophysical situations such as 
dark matter decay/annihilation, which will be studied in a separete paper. 
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Appendix A. Relevant Cross Sections 

In this appendix, we show the cross sections and energy loss rates adopted in this paper. 
Appendix A.l. Electron 

Appendix A. 1.1. Excitation of H We have adopted the almost same cross sections 
for atomic hydrogen and Coulomb losses as [1]. For the cross sections for the electron 
impact excitation of atomic hydrogen, we have adopted [7j. At high energy, we have 
used the Bethe approximation reviewed in [8], [9l [121 123] . 




cr2p(NR) 
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where Oq is the Bohr radius, Ry is the Rydberg energy and T = mv'^/2 represents the 
kinetic energy. The numerical coefficients are lnC2p = —0.89704, 72^ = 0.207985 and 
M|p = 0.55493. "NR" and "R" mean non-relativistic and relativistic respectively. 



a2.(NR) 



CT2s(R) 



A-koq 

T/Ry 

87ra; 



h2s + 



72s 



2s; 



where h 



2s 



mt^^/Ry 
0.11986 and 72, = -0.3125. 



(A.3) 
(A.4) 



T/Ry 



M|ln 



4C,T 



Ry 



0-n=3(R) 



^—Ml 



mf 2/Ry 



In 



(3' 



1-/52 



-/52 + lnC3 + 11.226^ 



(A.5) 
(A.6) 



where M| = 8.8989 x lO'^ and InCg = -0.2724. 

For the excitation to n > 3, we subtract Eq. flA.ip ~ Eq. flA.6p from total excitation 
cross section. 



aex(NR) 



0"ex(R) 



Anal 

T/Ry 



Mi In 





mw2/Ry 



In 



Ry y T/RyJ ' 



1-/32 



/3^ + lnC,^ + 11.2268 



(A.7) 
, (A.8) 



where InC^x = -0.5780,7ex. = -0.120575 and M^^ = 0.7166. 



Appendix A. 1.2. Ionization of H For the ionization cross section, we have adopted 
the following differential cross section. 



dai{E,e) 



A{E) 



de iHe/ey 0<e<-(E-/), (A.9) 

where E is the incident electron energy and e is the energy of the ejected electron. 
We choose e = 8 eV. In Eq. (lA.9p . a value of 2 is different from that of 2.1 originally 
suggested by [10]. This is why Eq. flA.9P can be analytically integrated. Two parameters 
A[E) and e are related to the total ionization cross section ai{E) 

<E) 



A{E) 
X{E) 



e 

E-I 



-[tan~^ X{E)]-\ 



(A. 10) 
(A.ll) 



where / is the ionization potential. The total ionization cross section for atomic hydrogen 
had been measured by |TT] in the range 14.6 — 4000 eV. Above 4000 eV, we used Bethe 
approximation [8], [12] . 



ai(NR) 



Airal 
Snan 



mf 2/Ry 



Ry 
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1-/32 



T/Rj 
-/3^ + lnCi + 11.2268 



(A.12) 
(A.13) 
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where = 0.2834, \nQ = 3.048 and 7^ = -1.6294 + ln(Ry/T). 

Appendix A. 1.3. e~ -H Collision For electron-hydrogen momentum transfer cross 
sections at low energies, we have adopted the results of [lH [Hj. The momentum loss 
cross section is described by 

^mt = 7^ $^[3 sin^(^r+i - Vr) + sm'(^ii - V^)] (A.14) 

where rj'^ and r]^ are the phase shift computed in fibi [16] . The cross section at 100, 
200 and 300 eV were calculated in [T7j. Cross sections at other energies were derived by 
interpolation and extrapolation. The energy loss due to electron-hydrogen momentum 
transfer is described by 

-dE] 2meE / . 

—J— = nHVe(Tmt[E), (A. 15) 

dt , m„ 

J rnt P 

where nip is the proton mass. 



Appendix A. 1.4- Coulomb Collision Incident electrons lose their energies due to elastic 
collisions with background electrons and photons. Energy loss is dominated by electrons 
at low energy since Coulomb cross sections are much larger than Compton cross sections. 
However, energy loss is dominated by photons at high energy since the number density 
of photons is much larger than that of electrons. 

For the energy loss due to Coulomb collisions with background electrons, we have 
adopted the following analytical formula [TSj . 

-dE' 



dt 



2.0X lO-V-97 / E-E, _i 



ci E^-^^ \E- 0.53E, 



eV-s-\ (A.16) 



where E is the incident electron energy in eV, E^ is the background electron energy in 
eV and Ue is the electron number density in cm~^. 

Appendix A. 1.5. Inverse Compton Scattering An important quantity which 
characterizes the behavior of inverse Compton scattering is jEcmb (photon energy in 
the electron's rest frame) where Eqmb is the energy of CMB photon. If '^Ecmb is much 
less than nie, Thomson scattering approximation is valid. Otherwise, Klein-Nishina 
cross section should be used. 

For inverse Compton scatterin with jEcmb mg, the energy spectrum of scattered 
photon is obtained by |1] in the limit /5 — > 1. However, we should keep (3 so as not to spoil 
the validity of Thomson approximation. After some tedious calculations, the number of 
collisions per unit time and unit scattered photon energy (E^) is given by 

d'^N 

—— = aTcn{EcMB)dEcMBf{x), (A.17) 
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where cry is the Thomson cross section and n{EcMB) is the differential number density 
of CMB photons and x = E^/Ecmb- The expressions for /(x) is given by 



167^/54 
2x In ( a; 
3 



1 + 



1 



X 



— 2x In 



1-/3 

1 1+/3 " 
xl -/3 



1 ^ 
for 



P 1-/3 
1-/3 



+ (x2(l + /3)-x(l 



1 + /3 
X 1 



< X < 1, 



1 + /3 1-/3 
for 1 < X < 



+ (x(l + /3)-x2(l 

1 + /3 
1-/3- 



(A. 18) 
(A.19) 



The coefficient is determined so that / f{x)dx is equal to unity. In the limit /3 — 1, 
Eq. ( ]A.19I) corresponds to the result of [1]. The number of collisions per unit time and 
the energy loss rate can be easily obtained. 

^ ^-dx = aTcn{EcMB)dEcMB = crrcncMB, 



dtdx 

d^N 4 

^^(E^ - EcMB)dx = -aTcn{EcMB)EcMBdEcMB7^P^ 

= ^cttcUcmbj'^P'^, 



(A.20) 



(A.21) 



where ucmb and Ucmb are the number and energy density of CMB photons. 

For inverse Compton scattering with ^Eqmb > f^ei "we should use Klein-Nishina 
cross section instead of Thomson cross section. The number of collisions per unit time 
and unit scattered photon energy is given by [22] 

2glng + (l + 2g)(l-g) + - ; , / , J l - g) 



dtdoi' ay 

X n{EcMB)dEcMB 



for a < a' < 



2 1 + 4a7g 
4a7^ 



(A.22) 



1 + 4a7 

where a = EcMB/nT'e, Oi' = E^/me and q = a' / ia'j'^ {1 — a' / 'j) . The number of collisions 
per unit time can be obtained by integrating Eq. (]A.22p . We shall assume that Eq. (]A.22p 
is valid for < g < 1, even though Eq. (]A.22p is quite invalid for < g < 1/47^. The 
contribution from the region < g < 1/47^ is 0(1/7^) and is negligible since Eqmb is 
much less than [22]. The number of collisions per unit time is given by [32] 



dtda' 



da' 



^ d^N 
dtdq 



dq = aTciJi{s)n{EcMB)dEcMB, 



where 



Ms) 



3 

2? 



9 



^ j ln(l + s) 



2s 



2 + 2s 



+ 4Li2(-s) 



and s = 4a'-f. The function Li2(x) is the dilogarithm 

ln(l - z) 



Li2(x) 



dz- 



(A.23) 



(A.24) 



(A.25) 
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The energy loss rate is given by 
f (fN 4 



(A.26) 



where 



Ms) 



^ + 6 + ^] Hl + s) 



6 + 13s + 8s2 + Ils3/i2 

OTsy 



+ 2Li2(-s) 



(A.27) 



Appendix A. 2. Photon 



Appendix A. 2.1. Photoionization The absorption of X-rays and 7-rays is studied in 
detail in [21]. Incident photons are mainly absorbed by hydrogen atoms and eject 
photoelectrons at low energies. The photoionization cross section for atomic hydrogen 
was reviewed in p!9| [20] . 

647rcrr f I \ exp(— 4?7 cot~^ 77) 
hu 



ai,(NR) 



exp{—27rri) 



(A.28) 



X 



4 

, 7(7- 



) [7^-lf/^ 



In 



7+ 



,7 



(A.29) 



7 + i \^ 27V/72 - 1 
where hu is the incident photon energy, a is the fine-structure constant, / is ionization 
energy, rj = Xj ^hv jl — \ and 7 = [hv + mf.)/me. The cross sections above are just 
halves of [121 [2D]- This is because there is only one electron in K shell in the case of 
hydrogen. 



Appendix A. 2. 2. Compton Scattering Incident photons interact with background 
electrons through Compton scattering. If photon energy is sufficient small, the energy 
of recoil electron is below the threshold energy of excitation and ionization of atomic 
hydrogen. Therefore the energy transferred to the recoil electron can be regarded as 
heating. Besides, a photon loses only a small fraction of its energy per scattering. The 
energy loss due to Compton scattering is described by [2T] 
^-dE' 



dt 



meUeCaTX g[x) 



(A.30) 



Compton 



where x = hu/me and 



[X 



3)(x+ 1 



+ 



2(3 



x^ 
17x- 



ln(l + 2x) 
31a;2 + 17x3 - ■ 



3 



x" 



2x)' 



ix^ 



for X ^ 1: 



In 2x 

6 



for X ^ 1; 



(A.31) 
(A.32) 
(A.33) 
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If photon energy is as large as the electron mass, a photon loses a sizable fraction of its 
energy per scattering. In this case, it is necessary to calculate the energy distribution of 
the recoil electrons. The cross section is given by the following Klein-Nishina formula: 



da 
de 



{hv) 



3(Tt rUe 

8 (H- 



hu 

e 



e 

hu 



e 



me 
hu 



rrie 
e 



hu 



for 



nip 



nif, + 2hu 

where e is the scattered photon energy. 



-hu < € < hu. 



(A.34) 



Appendix A. 2. 3. Pair Creation If photon energy is larger than 2me, it is possible 
to create electron-positron pair. The energy and momentum conservation, however, are 
only possible if another particle is present. The differential cross section for pair creation 
in nuclei is given by the Bethe-Heitler formula [IHl [21] : 



dE. 



ar^ 



-- - 2E+E J ^ • ; + m 
P+P 



+ L 



[ElEl+pW 



o%p'i 
8 E^E^ 



p-1 



mlEj 



pI 



X 



P%pt 

E.E^ 



3 p+p- '^P+P- 



p-L^ , E+E_-pI^ 2E^E+E_ 



P'l 



P+ 



where 



P± 



U 



El. - ml, 



In 

-- In 



E+E_ + p+P- + TTlg 

E.^.E^ — p+p- + ml 
E±+p± 



P+P- 



(A.35) 

(A.36) 
(A.37) 

(A.38) 



E± -p± 

and E± is the energy of positron (electron) energy. The analytical expression for cross 
section is given by [211 [2S] 

a 



1 I 28 
m' < — In 2A; 
1 9 



218 / 2 



7 



61n2A; - - 
2 



2 , 
- In^ 2k 
3 



\xi^2k 



7r 



ln2A; 



TT 

29 



2C(3) 



9 X 256 



ln2A; 



2 

k 

77 



27 X 512 



16 

+ 



ln2A; 



for k> A 



(A.39) 



where k = E^/m^. Convenient approximate formulas are given by [26] which are valid 
for k < 20, 



^ = ar20o^[l + 0.135( 
ax 



0.52)^(1 - z% 



where 



X 



E+ - me 
E.y — 2mp 



(A.40) 
(A.41) 
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z = 2^/x{l~x), (A.42) 

and 00 is the differential cross section for equal partition of energy, = E^ = E^/2. 
The second term in the square bracket should be dropped when it becomes negative 
(below k = 4.2). (po is given by 



;i-7i) 



-(4 - 7i)(Li - 1) - 7iai("i - 1) - 7i"i(-^i - "i) 



where 



7i 



2 



1 



In 



k 

2 + 



We get from Eq. ( 1A.40I) for the total cross section 

^ 2 

—or 
4 

ar,2(O.7760o + 0.0180^) for k > 4.2. 



for k < 4.2, 



(A.43) 

(A.44) 
(A.45) 

(A.46) 

(A.47) 
(A.48) 



Appendix A. 2. 4- Photon-Photon Scattering If the photon energy is below the effective 
threshold energy of the double photon pair creation, photon-photon scattering (77 — > 
77) process becomes significant. The photon-photon scattering rate for E^Ecmb < ''^e 
is given by [27] 

6 / 77. \ 3 



P{E, 



3.33 X 10 



11 



-CMB 



Normalized distribution of secondary photons of energy E' p{E' E^), is given by 



(A.49) 



20 1 



E^ 



E^ 



for < < E^. 



The distribution p{E'E^) satisfies 



p{E' E^)dE' 



p{E'E,)E'dE'=E,. 



(A.50) 



(A.51) 



(A.52) 



The above formulas are not valid for a larger value of E^. However, photon-photon 
scattering is not significant for high energy photons since double photon pair creation 
is the dominant process. Therefore, instead of using exact formulas, we simply neglect 
photon-photon scattering for E^yEcuB > f^l- 



Electron and Photon Energy Deposition in Universe 



20 



Appendix A. 2. 5. Double Photon Pair Creation For high energy photon, double photon 
pair creation (77 —>■ e^e~) is the dominant process. The total cross section for double 
photon pair creation is given by [28] 

1 



a 



(3-/3^)lnl±|-2/3(2-/?2) 



(A.53) 



where /3 is the electron (positron) velocity in the center-of-mass system. The relationship 
between (3 and E^, Ecmb and 9 which is the angle between the momenta of the colliding 
photons is easily obtained: 



E^Ei 



y^CMB 

2ml 



cos I 



(A.54) 
(A.55) 



Clearly, the threshold energy for double photon pair production is E^ = ml/ Ecmb, 
head-on photon collision = tt, s = 1). For calculation of the absorption probability, we 
should average the above cross section over the distributions for isotropically distributed 
photons [28], 



l-l~2ml/E^EcMB 
2/1 
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COS 6*) ad COS ^ 
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\yiuj - 2 



y 



-(l + 2tO 



where 



E^Ei 



•y^CMB 



1 > 0, 



Differential spectra of electrons and positrons are given by [29] 



(A.56) 

(A.57) 
(A.58) 



da 
IK 



1 4 

■nr^m^ 

^^^^CMB 



AE"^ AEcMBEeEp 8E, 
m ■ 



CMB 



E 



EeEp 



miE 



mz 



^ 2{2EcMBE-ml)E^ ^ ^ 
mlEeEp 



mz 



EcmbE J E^E^ 
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where Ee{Ep) is the energy of electron (positron) and E is the total ener; 
Ee + Ep = E^ + Ecmb- The limits of the variation of Ee{Ep) is given by 




(A.59) 

lY, E = 

(A.60) 
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Appendix A. 3. Positronium 



Here, we show the energy spectrum from three-photon annihilation of positronium. The 
energy spectrum is continuous, as allowed by conservation of momentum. It has been 
calculated in [31] with photon energy rj normalized by electron mass. 
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2(1 -r?) 
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ln(l — rj) 
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The function F{ri) is normalized so that 



drjF^f]) = 1. 
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(A.61) 



(A.62) 
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